The Neutrino Response of Low-Density Neutron Matter from the Virial Expansion 
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We generalize our virial approach to study spin-polarized neutron matter and the consistent 
neutrino response at low densities. In the long-wavelength limit, the virial expansion makes model- 
independent predictions for the density and spin response, based only on nucleon-nucleon scattering 
data. Our results for the neutrino response provide constraints for random-phase approximation or 
other model calculations, and we compare the virial vector and axial response to response functions 
used in supernova simulations. The virial expansion is suitable to describe matter near the supernova 
neutrinosphere, and this work extends the virial equation of state to predict neutrino interactions 
in neutron matter. 
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I. INTRODUCTION 

Neutrinos radiate 99% of the energy in core-collapse 
supernovae. The scattering of neutrinos and the physics 
of the explosion are most sensitive to the properties of 
low-density nucleonic matter P, Q , which is a complex 
problem due to strong coupling with large scattering 
lengths, clustering in nuclear matter and the non-central 
nature of nuclear interactions. For low densities and high 
temperatures, the virial expansion provides a tractable 
approach to strong interactions, and in previous works we 
have presented the virial equation of state of low-density 
nucleonic matter '3, The predicted large symmetry 
energy at low densities has been confirmed in near Fermi 
energy heavy-ion collisions 0. 

The virial approach can be used to describe mat- 
ter in thermal equilibrium around the neutrinosphere 
in supernovae. The temperature of the neutrinosphere 
is roughly T ~ 4MeV from about 20 neutrinos de- 
tected in SN1987a 0, 0, and the density follows from 
known cross sections of neutrinos with these energies 
n ~ 10^"'^ — lO-'^^g/cm^. For neutron matter, the virial 
expansion in terms of the fugacity z = e^^'^ is valid for 



A3 



^ + 0{z^) < 4 • 10" (T/MeV)3/2 g/cm^ , (1) 



where we require z < 1/2 and A denotes the thermal 
wavelength A = (27r/rnT)^/^. Therefore, the virial ap- 
proach makes model-independent predictions for the con- 
ditions of the neutrinosphere, based only on the experi- 
mental scattering data. 

In this paper, we use the virial expansion to describe 
how neutrinos interact with low-density neutron mat- 
ter. We focus on neutral-current interactions, and leave 
charged-current reactions and nuclear matter to future 



works. Our long-term goal is a reliable equation of state 
and consistent neutrino response for supernovae. 

The free cross section per particle for neutrino-neutron 
elastic scattering is given by 



1 rfcTQ 
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+ C,2(l + cos0) , (2) 



where Gp is the Fermi coupling constant, E,y the neu- 
trino energy, and 6 the scattering angle. The weak axial 
coupling is Ca = ga/'2, with ga = 1.26 the axial charge of 
the nucleon. The weak vector charge is Ci, = — 1/2 for 
scattering from a neutron. Eq. ||2Jl neglects corrections of 
order E^/m from weak magnetism and other effects 

In the medium, this cross section is modified by the 
vector response S^ici) and the axial response S'a('?) 
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where Sv and Sa describe the response of the system 
to density and spin fluctuations respectively, and q = 
2E„ sm{6/2) denotes the momentum transfer. We will 
discuss the approximations for Eq. ^ in Sect. Ill CI In 
the following, we will use the virial expansion to provide 
model-independent results for the response in the long- 
wavelength {q — !■ 0) or forward-scattering limit. 

This paper is organized as follows. We extend the virial 
equation of state to spin-polarized matter in Section ^ 
and derive the consistent long-wavelength response. Fur- 
ther details on the virial equation of state can be found 
m Rcfs. Hi. In Section IIIII we present results for the 
spin virial coefficients, the pressure and entropy of spin- 
polarized neutron matter, and the neutrino response. We 
compare our results to Brueckner calculations, and to 
random-phase approximation (RPA) response functions. 
Finally, we conclude in Section Hvl 
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II. FORMALISM 

The virial expansion is a general, model-independent 
approach for a dilute gas, provided the fugacity is small 
and for temperatures above any phase transitions. Under 
these conditions, the grand-canonical partition function 
can be expanded in powers of the fugacity. The second 
virial coefficient 62 describes the term in this expan- 
sion and is dire ctly related to the two-body scattering 
phase shifts [lOl . The relation of the third virial co- 
efficient to three-body scattering is not straightforward, 
and was only studied for special cases |0, 0, 0| . The 
virial expansion is not a perturbative k-pa^ expansion, 
and its great advantage is that it includes bound states 
and scattering resonances on an equal footing. 



as a function of the spin- up one A,T), and gen- 

erate the virial equation of state in tabular form for a 
range of z+ values. This maintains the thermodynamic 
consistency of the virial equation of state. 

Finally, we will also discuss results for the entropy. The 
entropy density s = SjV follows from differentiating the 
pressure with respect to the temperature s = (pTP^iii ■ 
This leads to 

5^ 1 1 

s = — -n+logz+-n_ logz_ 

+ J(fe:^,l(4 + ^-)+2 6;,0^+2-): (8) 

where 6'(T) = db{T)/dT denotes the temperature deriva- 
tive of the virial coefhcients. 



A. Spin-Polarized Matter 

The virial equation of state is easily generalized to spin- 
asymmetric systems. For two spin components, we de- 
note the chemical potential for spin up and spin down 
particles by and with fugacity z+ — e^+l^ and 
z_ = e^-l'^ respectively. For the virial equation of state 
we expand the pressure in a power series of the fugacities 

T 

P = (z+ + 2- + 6n,i (z+ + ) -I- 2 6„,o z+z- + 0{z^)) . 

(4) 

The second virial coefficients for like spins and 6„^o 
for opposite spins are related to the two-particle partition 
function and are given in terms of the scattering phase 
shifts in the next section. The densities follow from dif- 
ferentiating the pressure with respect to the fugacities. 
For the density of spin-up neutrons = {dp,_^P)T = 
Z+/T [dz^P)T we thus have 

n+ = ^ (z+ + 2 6„,i zl + 2 6„,o z+z^ + 0(2^)) , (5) 

and likewise for the density n_ of spin-down neutrons 

n_ = ^ (z_ + 2 5„.i zl+2 bn,o z-z+ + 0{z^)) . (6) 

The total density n and the spin polarization A are then 
given by 

n = n+ + n_ and A = . (7) 

n+ + n_ 

In this work, we truncate the virial expansion after sec- 
ond order in the fugacities. This leads to an equation of 
state that is thermodynamically consistent. 

The dependence of the total density and the spin 
polarization on z^ and 2_ can be inverted to yield 
the virial equation of state directly in terms of 
P(z+(n, A, T), z_(n. A, T), r). In practice, for a given 
spin polarization, we determine the spin-down fugacity 



B. Spin Virial Coefficients 

The second virial coefficient 6„.i describes the interac- 
tion of two neutrons with the same spin projection. To 
this end, we generalize the second virial coefficient of the 
spin-symmetric system 0, 0, 0| to 

9I/2 /-oo 

6„,i(T) = ^/ dEe-^/'^6T\E)~2-'/\ (9) 
ttT Jo 

where — 2~^/^ is the free Fermi gas contribution and 
6\°^{E) is the sum of the isospin and spin-triplet elas- 
tic scattering phase shifts at laboratory energy E. This 
sum is over all partial waves with angular momentum L 
and total angular momentum J allowed by spin statistics, 
and includes a degeneracy factor {2J + 1)/(2S' + 1), 

5T\E)^Y.^-^^-^^AE) 

L,.J 

= i (53p„ -f ,53p^ + ^ J3p^ . . . (10) 

The factor 1/(25' -I- 1) = 1/3 arises because the same 
spin projection, e.g., for up spins Ms = -1-1, is 1/3 of the 
possibilities Ms — —1,0,1. Note that we have neglected 
the effects of the mixing parameters due to the tensor 
force. We expect that their contributions are small for 
low densities. 

Two neutrons with opposite spin projections have a 
probability 1/2 to be in spin 5 = or 5* = 1 states, 
thus the second virial coefficient for opposite spins 6„.o is 
given by 

9I/2 fCO 

bnAT)^^ I dEe-^/^^Sr{E), (11) 

where 6q°^{E) is the sum of allowed isospin-triplet elas- 
tic scattering phase shifts with degeneracy factor (2 J -t- 
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l)/(2(2^+l)), 

= 2 ^ 6 ° ^ 2 ^ 6 ' ^ 2 °' ^ 



(12) 



The second virial coefficient for spin-symmetric neu- 
tron matter 6„ is the sum over hke and opposite spins, 



(13) 



and consequently the sum of the total phase shifts 
given above determines 6„ with d*°^{E)/2 = Sq°^{E) + 
S\°^{E) [23. In addition, we define the axial spin virial 
coefficient 6q as 

ba = bn,i - &„,o ■ (14) 
Thus, if only S-wave interactions are present, one has 

bn,i = -2-'/' and ba = -bn - 2-3/2 ^^^^ 

C. Neutrino Response 

Neutrino scattering from a many-body system can 
be expressed in terms of the vector Sv{q,w) and axial 
Saigjiv) dynamical response functions. These describe 
the probability for a neutrino to transfer momentum q 
and energy w to the medium. Integrating over energy 
transfer, we define the static vector S.v{q) and axial Sa{q) 
response functions 



dwSy^aiq.w) 



(16) 



Here scattering kinematics limits the energy transfer to 
be space-like \w\ < q. At low densities nucleons are non- 
relativistic, and therefore we expect the vector response 
to have little strength in the time-like region so that 



Sviq) 



dw Sv (q, w) . 



(17) 



The axial response can have contributions from multi- 
pair states in the time-like region even in the long- 
wavelength limit due to non-central nuclear interac- 
tions [l5j. However, neutron matter at very low density 
can be described using a pion-less effective field theory 
where non-central interactions are sub-leading. There- 
fore, in this paper we approximate the axial response by 



The static structure factor for the density response is 
then given by 



3j. giq r 



where the sum is over all many-body eigenstates \ j) with 
energy Ej, the partition function is ^ = ^je~^^^ and 
(3 — 1/T. For the spin response, the density operator is 
replaced by the spin density tp^ {r)(Tip{r). 

In the long- wavelength limit, the vector response of 
the spin-symmetric system is given by (see also Ap- 
pendix B in 



S4q = 0) = 



T 



[dP/dnYj 



(20) 



For the symmetric system, the total chemical potential 
is /i = (/i+ -I- /i_)/2, with fugacity z = ^fz+zZ = e^l^ ^ 
and the virial equation of state (see also Ref. 0) yields 
for the consistent vector response. 



n \ oz 



1 + AbnZ 
1 -I- 2b„z 



(21) 



Following Burrows and Sawyer |17| , we define the spin- 
difference or axial chemical potential fj,a = itJ-+ — M-)/2 
and the axial fugacity Za = yj zj^j z-. The axial response 
of the spin-symmetric system is then given by 



5a(g = 0) 



d 



("4 



n dza 

and the virial expansion, Eq. I^J, leads to 

2baz 



z„=l 



Sa{q = 0) = 1 



1 + 2bnZ 



(22) 



(23) 



The long-wavelength limit of the axial response is also 
related to the spin susceptibility x. 



Sa{q = 0) 



XF 



nT 



ln.T.A=0 



(24) 



where / denotes the free energy density, and xf = Mn'^/^^ 
is the spin susceptibility of a free neutron gas, with the 
neutron magnetic moment Finally, the response 

functions are normalized to unity in the low-density limit 
S'„(0) = 5a(0) = 1 for z = n = 0. 



III. RESULTS 



A. Spin Virial Coefficients 



Saiq) 



dw Sa{q, w) . 



(18) 



We first calculate the virial coefficients 6„ and ba from 
the T = 1 np phase shifts obtained from the Nijmegen 



TABLE I: The second virial coefficient 6„ and the axial virial 
coefficient ba for different temperatures. The results labeled 
CIB take into account the effects due to charge-independence 
breaking (CIB) on the scattering length with a„„ — —18.5 fm. 
We estimated an error of < 5% for the higher temperatures 
T ^ 25 MeV due to the truncation of the integration over the 
phase shifts at E ^ 350 MeV. 



T[MeV] 


bn 


with CIB 




ba 


with CIB 




I.OO 


0.288 


0.251 


0.032 


-0.641 


-0.604 


-0.031 


2.00 


0.303 


0.273 


0.012 


-0.655 


-0.625 


-0.007 


3.00 


0.306 


0.279 


0.004 


-0.655 


-0.629 


0.006 


4.00 


0.306 


0.283 


0.001 


-0.652 


-0.628 


0.014 


5.00 


0.306 


0.285 


0.000 


-0.648 


-0.627 


0.020 


6.00 


0.306 


0.286 


0.001 


-0.644 


-0.624 


0.023 


7.00 


0.307 


0.288 


0.002 


-0.640 


-0.621 


0.026 


8.00 


0.307 


0.289 


0.004 


-0.637 


-0.619 


0.028 


9.00 


0.308 


0.291 


0.007 


-0.634 


-0.616 


0.029 


10.00 


0.309 


0.292 


0.009 


-0.631 


-0.614 


0.029 


12.00 


0.310 


0.295 


0.013 


-0.625 


-0.610 


0.029 


14.00 


0.313 


0.299 


0.017 


-0.621 


-0.607 


0.028 


16.00 


0.315 


0.302 


0.020 


-0.617 


-0.604 


0.026 


18.00 


0.318 


0.305 


0.022 


-0.614 


-0.602 


0.024 


20.00 


0.320 


0.308 


0.023 


-0.612 


-0.600 


0.021 


22.00 


0.322 


0.311 


0.023 


-0.610 


-0.598 


0.019 


24.00 


0.324 


0.313 


0.022 


-0.608 


-0.597 


0.018 


26.00 


0.326 


0.315 


0.021 


-0.607 


-0.596 


0.017 


28.00 


0.327 


0.317 


0.018 


-0.606 


-0.595 


0.016 


30.00 


0.329 


0.318 


0.015 


-0.605 


-0.595 


0.016 


35.00 


0.330 


0.321 


0.004 


-0.602 


-0.593 


0.018 


40.00 


0.330 


0.321 


-0.009 


-0.599 


-0.591 


0.024 


45.00 


0.328 


0.319 


-0.025 


-0.596 


-0.588 


0.031 


50.00 


0.324 


0.316 


-0.041 


-0.592 


-0.584 


0.041 



partial wave analysis [iSj . This neglects the small charge 
dependences in nuclear interactions. We have included all 
partial waves with L ^ 6. For the higher temperatures, 
T ^ 25 MeV, there is a < 5% error due to the truncation 
of the integration over the phase shifts at E ^ 350 MeV 
(the extent of the partial wave analysis). This error was 
estimated by assuming constant total phase shifts and 
varying the energy cutoff to E > 350 MeV. 

Our results for the virial coefficients and their tempera- 
ture derivatives Tb'{T) are listed in Table|l| As discussed 
in Ref. 4], the virial coefficients are dominated by the 
large S-wave scattering length physics (a„p — —23.768 fm 
and ann = — 18.5fm), but effective range and higher par- 
tial wave contributions are noticeable. For example, for 
T — 5 MeV, the virial coefficients obtained only from a„p 
arc bn = 0.44 |1| and ba = -bn - 2-^/2 = -0.80. In the 
unitary limit where the scattering length ja^j = ±oo and 
S{E) = 7r/2, the second virial coefficients are independent 
of the temperature and given by bn = 3/2^/2 = 0.53 [l| 
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FIG. 1: (Color online) The pressure P versus density n for 
T = 20 MeV and various spin polarizations A = (symmet- 
ric), 0.5 and 1 (fully polarized). We also compare our results 
to Brueckner calculations of Bombaci et al. (crosses with dot- 
ted lines) |23] for A = 1,0.5,0 (top to bottom). The circle 
indicates where the fugacity is z = 0.5 for A = 0, and for the 
other spin polarizations the virial curves end at 2: = 1. 

and ba — —5/2^/^ = —0.88. Therefore, the virial expan- 
sion is well defined for resonant interactions, in contrast 
to the kpas expansion. 

We find that the second virial coefficients are approx- 
imately independent of temperature over a wide range, 
and consequently Tb'{T) « 0. As a result, the ther- 
modynamic properties of spin-polarized neutron matter 
and the long-wavelength response scale as a function of 
the fugacities, which depend on density and tempera- 
ture through Zi(n+/r^/2^^_y'7^3/2-j j^j. ^ _ _|_ ^^j^^ i = —. 

This scaling can also be expressed in terms of the Fermi 

2/3 

temperatures Tf.i ~ i sjid thus the properties of neu- 
tron matter scale with T/Tp^, ~ T/7\/^ only. In Ref. [3| 
we found that spin-symmetric neutron matter scales to 
a very good approximation. The virial scaling symme- 
try is exact for cold atomic gases tuned to a Feshbach 
resonance pol | and has been verified experimentally by 
Thomas et al. 

In Table we also study the effects of charge- 
independence breaking (CIB) on the scattering length. 
We estimate CIB effects as discussed in Ref. 0. CIB for 
the virial coefficients is largest for T < 5 MeV and leads 
to a 10% reduction in magnitude of the virial coefficients. 



B. Pressure and Entropy of Spin- Polarized Matter 

We have previously found 0] good agreement of the 
virial equation of state for spin-symmetric neutron mat- 
ter with microscopic Fermi hypcr-ncttcd chain (FHNC) 
calculations of Friedman and Pandharipande 22] for den- 
sities up to n < rio/lO, where no = 0.16 fm^^ is the sat- 
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FIG. 2: (Color online) The entropy per particle S/N versus 
density n for T — 20 MeV and various spin polarizations A = 
(symmetric), 0.5 and 1 (fully polarized). We also compare 
our results to Brueckner calculations of Bombaci et al. (crosses 
with dotted lines) Q for A = 0, 0.5, 1 (top to bottom). The 
circle indicates where the fugacity is z = 0.5 for A = 0, and 
for the other spin polarizations the virial curves end at z — 1. 



uration density of symmetric nuclear matter, and pub- 
lished temperatures T ^ 10 MeV. For nuclear matter, 
the FHNC results fail to describe clustering with alpha 
particles at low densities . 

Our virial results for the pressure and entropy of spin- 
polarized neutron matter are shown in Figs. ^ and |21 for 
T — 20 MeV and polarizations A = (symmetric), 0.5 
and 1 (fully polarized). For this temperature, we can 
compare the virial results to Brueckner calculations of 
Bombaci et al. As shown in Fig. [3 the Brueckner 

entropy agrees well with the virial results. For the pres- 
sure, the effects of a spin polarization are smaller, and 
in addition there is some uncertainty in the Brueckner 
calculations, since the pressure was obtained from the 
energy by a numerical derivative. Before we discuss the 
neutrino response, we note that it is difficult to calcu- 
late the long-wavelength response at low densities from 
the Brueckner or FHNC results, since the response is ob- 
tained by differentiating the pressure. 



C. Neutrino Response 
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FIG. 3: (Color online) The vector and axial response of 
neutron matter for T — 4 MeV. In addition to the long- 
wavelength virial response, we also show the RPA response of 
Burrows and Sawyer 17] for neutron matter and various mo- 
mentum transfers q — 0, 3T and 6T. For this density range, 
the fugacity in the virial expansion is z < 0.42. 



lowest order in the density. With z « nX^ /2, we have 



> 1. 



(25) 



since bn = 0.31 from Table In a Landau- Fermi liquid, 
the vector response is given by S'„(0) = 1/(1 + Fq) > 1 
for neutron matter, where the Landau parameter for the 
density-density interaction is i^o < ||24| . 

In contrast, the spin-spin interaction is repulsive (this 
follows from the Pauli principle, when the density-density 
interaction is attractive), and the virial axial response 
gives S'a < 1 for low-momentum transfers. This is seen 
in the low-density limit. 



5a(g-0)«l + 6anA3<l, 



(26) 



where ba — —0.65 from Table Analogous to the vec- 
tor response, the axial response for a Landau-Fermi liq- 
uid is given by S'a(O) = 1/(1 + Go) < 1 for neutron 
matter, since the Landau parameter for the spin-spin 
interaction is Go > |2^. Although the virial densi- 
ties and temperatures are not in a Fermi liquid regime 
{z ~ (Tf/T)^/2 < 1), the deviation of the vector and 
axial response from a free gas is determined by nuclear 
interactions, and thus is the same for low and high tem- 
peratures. 



Our virial results for the long-wavelength vector and 
axial response are presented in Fig. |31 The neutron- 
neutron interaction is attractive at long distances and 
thus increases the probability to find two neutrons close 
together compared to a free neutron gas. These density 
fluctuations increase the local weak charge and produce 
a vector response Sy > 1 for low-momentum transfers. 
This is easily seen by expanding the vector response to 



D. Comparison to RPA calculations 

Most present calculations of the neutrino response are 
based on the random-phase approximation (RPA) 0, 
ITtI I25I I , which gives the linear response of a mean- field 
ground state to neutrinos. The RPA response thus ne- 
glects clustering and is incorrect for nuclear matter at 



6 




FIG. 4: (Color online) The long-wavelength virial response ^1^. 5: (Color online) The total response of neutron matter 
of neutron matter for T = 4 MeV is compared to the RPA gi^en by the ratio of the total cross section for elastic neutrino- 
response with effective-range expansion (ERE) contact inter- neutron scattering in the medium compared to free space, 
actions and average neutron energy (E) = 3T/2. The results shown are for T = 4 MeV and neglect the small 

momentum dependence of the vector and axial response. 



subnuclear densities |26l |. Since there is no clustering 
in neutron matter, a comparison of the virial with RPA 
response assesses the interactions used in present RPA 
calculations, as well as the random-phase many-body ap- 
proximation for low densities and high temperatures. 

As an example, we compare our virial results to 
the nonrelativistic RPA calculations of Burrows and 
Sawyer where the RPA interaction is chosen to re- 
produce Landau-Fermi liquid parameters for symmetric 
nuclear matter. We compare to the approach of Ref. ^3 j 
because these results have been used in supernova sim- 
ulations and they are somewhat simpler and thus 
more transparent than Refs. I2,'t|| . Since Burrows and 
Sawyer do not present results for pure neutron matter, 
we have calculated the RPA response following Ref. 0| . 
For completeness, we give the necessary equations in Ap- 
pendix^ For low-density neutron matter, the effective 
mass is well approximated by the free mass and we 
thus use m* /m = 1. 

In Fig. 131 we compare the RPA results for T = 4 MeV 
to the virial response. We find that the RPA axial re- 
sponse is repulsive {Sa < 1) and on a qualitative level 
similar to the virial response. However, the RPA vector 
response is also repulsive, in contrast to our virial result. 
This is because Ref. uses Landau parameters of sym- 
metric nuclear matter for all proton fractions. In partic- 
ular. Burrows and Sawyer use for the spin-independent 
part of the interaction Fq -|- Fg Ti • T2, with Fq = —0.28 
and Fq = 0.95 0|, and the density-density interaction 
for neutrons is thus Fo{T = 1) = Fq + F^ = 0.67. This 
makes the incorrect assumption that induced interactions 
in nuclear and neutron matter are identical. For the virial 
coefficient 6„, the total phase shift is attractive 0. This 
leads to an attractive vector response {Sy > 1) at low 
densities and low- momentum transfers. The RPA results 
of Refs. 0, 123 have an attractive vector response, con- 



sistent with the mean-field equation of state. However, 
the axial interaction of Ref. [23 is not constrained at the 
mean-field level and may be more poorly determined. 

The RPA provides a model to study the momentum 
dependence of the response functions. For a neutrino 
with energy Ei, = ZT, the maximum momentum transfer 
is (Zmax = 2i?^ = 6T. In addition to the long-wavelength 
response, FigOshows the RPA results for various momen- 
tum transfers. This demonstrates that the RPA response 
has a very weak momentum dependence. Consequently, 
the long- wavelength response provides strong constraints 
for all relevant momentum transfers. 

Next, we calculate the RPA response, when we use con- 
tact interactions that are constrained by nucleon-nucleon 
scattering. In order to obtain cutoff-independent re- 
sults and correctly include the large scattering length 
and effective range at low density, it is necessary to sum 
particle-particle ladders and work with the T matrix (see 
also 13). This leads to Landau parameters /o+ffo cri •o'2 
(-Fb = rnkF/o/ir'^ and Gq = n^kpgo/ir'^) for the antisym- 
metrized interaction with 



/o = 



27r/m 



l/a„p - mre E/2 



and go = - fo-. 



(27) 



where = 2.68 fm is the effective range and E de- 
notes the relative energy. In order to make a sim- 
ple estimate, we take an average relative energy (E) = 
{{Pi ~P2f') /4:ni — 3T/2 and calculate the RPA response 
with these Landau parameters. The resulting vector and 
axial response is shown in Fig. ^ The axial response 
agrees nicely with our virial result, but for the vector re- 
sponse there is only a good agreement at low densities. 
The differences at higher densities could be due to us- 
ing an average energy, since the vector response is more 
sensitive to the latter. 
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Finally, we show the total response of neutron mat- 
ter for T — 4MeV in Fig. O The total response is 
given by the ratio of the total cross section for elas- 
tic neutrino-neutron scattering in the medium compared 
to free space. We neglect the small momentum depen- 
dence of the vector and axial response, and thus have 
ct/^o = {&9lSa{0) + 25„(0))/(6g2 + 2). We find for ex- 
ample a factor 0.72 reduction of the total response at 
n = 0.0016 fm^'^ rio/lOO. This is 10% larger compared 
to the RPA response of Burrows and Sawyer. 



IV. CONCLUSIONS 

We have extended our virial approach to study 
spin-polarized neutron matter and the consistent long- 
wavelength response. The virial expansion is suitable 
to describe matter near the supernova neutrinosphere, 
and this work extends the virial equation of state 0, Q 
to predict neutrino interactions in neutron matter. Our 
results include the physics of the large neutron-neutron 
scattering length in a tractable way. We have found that 
the spin virial coefficients are approximately temperature 
independent over a wide range. The properties of spin- 
polarized neutron matter and the response therefore scale 
with density and temperature as discussed in Ref. .4] • 

The virial expansion was used to make model- 
independent predictions for the pressure and entropy of 
spin-polarized matter, for the vector and axial response, 
and the cross section for neutrino-neutron scattering in 
the medium. The virial pressure and entropy of spin- 
polarized neutron matter are similar to Brueckner re- 
sults 23] , but the virial approach has a well-defined range 
of validity and is directly based on scattering data. 

The virial equation of state predicts an attractive vec- 
tor and a repulsive axial response in the long-wavelength 
limit at low densities. The total neutrino response is 
suppressed in matter compared to the response of a free 
neutron gas. This provides a benchmark for many-body 
calculations of the response functions. As an example, 
our results for the neutrino response disagree with the 
RPA response of Burrows and Sawyer due to the in- 
teraction model used for the latter. The RPA was used to 
study the momentum dependence of the response func- 
tions. We have found a very weak dependence on mo- 
mentum transfer (independent of sign and magnitude of 
the interaction). We therefore conclude that the long- 
wavelength virial response provides strong constraints at 
low densities for all relevant momentum transfers. 

Important areas of future work are the extension of 
these techniques to charged-current interactions and to 
the neutrino response in nuclear matter |26l |. In ad- 
dition, a generalization of the virial expansion beyond 



Section III CI may offer insights to the effects of multi- 
pair states on the long-wavelength response at low den- 
sities The third virial coefficient can be used to 
provide error estimates For the neutrino response, 
an effective field theory calculation of the dominant large 
scattering length contributions to the third spin virial co- 
efficients would be very useful. 
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APPENDIX A: RPA RESPONSE 

The static structure function S{q) is given in terms of 
the polarization function ^) by 



^il) = - f dtd^ 

TT J 1 



where the polarization function in RPA reads 

n°(<z,^) 



(Al) 



(A2) 



The Landau interaction used by Burrows and Sawyer [l7| 
is vo = fo ~ 1.76 • lO^^MeV^^ for the density response 
and vq = go — 4.50 • 10~^ MeV~^ for the spin response. 
The real and imaginary parts of the free polarization 
n''((7,tj) are derived in the Appendix of Ref. [T^ . 



Ren°(g,w) = 



27r2g/3 



ds 

— In 

s 



1 



1 + e^^t^-is-Qy 



Imn°(g,tj) = 



2TTq(3 



In 



1 + 6' 



(A3) 
(A4) 



with Q — y/nif3/2 {—lo / q + q/ {2m)). Finally, the density 
is given by 



(27r)3 1 -He'3(pV(2"i)-/.) 



(A5) 



which determines the chemical potential for 
Eqs. IIA3IA4II . 
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